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Abstract 



Analytical expressions for the width and conductance peak distributions of 
irregularly shaped quantum dots in the Coulomb blockade regime are pre- 
sented in the limits of conserved and broken time-reversal symmetry. The 
results are obtained using random matrix theory and are valid in general for 
any number of non-equivalent and correlated channels, assuming that the un- 
derlying classical dynamic of the electrons in the dot is chaotic or that the 
dot is weakly disordered. The results are expressed in terms of the chan- 
nel correlation matrix which for chaotic systems is given in closed form for 
both point-like contacts and extended leads. We study the dependence of the 
distributions on the number of channels and their correlations. The theoreti- 
cal distributions are in good agreement with those computed in a dynamical 
model of a chaotic billiard. 

PACS numbers: 73.40. GK, 05.45+b, 73.20.Dx,24.60.-k 
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I. INTRODUCTION 



One of the most interesting aspects of electron transport in submicron scale devices 
is the interplay between quantum coherence and aperiodic but reproducible conductance 
fluctuations. Over the past decade the phenomenon of universal conductance fluctuations in 
disordered systems (where impurity scattering dominates) has been understood through the 
use of stochastic models. More recently, a new generation of experiments^ was designed to 
measure conductance fluctuations in the ballistic regime where the dynamics of the electrons 
in the device is determined by the geometry of its boundary. The stochastic approach to 
these systems is justified by the underlying classical chaotic dynamics. This situation is 
distinct from the diffusive case, where the corresponding classical limit of the quantum 
problem is not fully understood. 

In this paper we discuss the conductance fluctuations in quantum dots.H These are 
semiconductor devices in which the electrons are confined to a two-dimensional region whose 
typical linear dimension is in the submicron range. In particular we are interested in the 
Coulomb blockade regime where the leads are weakly coupled to the dot, either because the 
leads are very narrow, or due to the presence of potential barriers at the lead-dot interface.! 
The electrons inside the dot are characterized by isolated resonances whose width is smaller 
than their average spacing, and conductance occurs through resonant tunneling. As a con- 
sequence, the conductance peaks when the Fermi energy matches a resonance energy of 
the electrons inside the dot and an additional electron tunnels into the dot. Such a sys- 
tem resembles the compound nucleus in its region of isolated resonances.! The macroscopic 
charging energy required to add an electron to a dot is determined by its capacitance C and 
is given by e 2 /C. Since C is a constant which is determined essentially by the geometry 
of the dot, the conductance exhibits equally spaced oscillations as a function of the gate 
voltage (or Fermi energy). At low temperatures T <C kT < A the width of the conductance 
peaks is ~ kT, but the heights exhibit order of magnitude variations.! - ! 

When the electron-impurity mean free path is larger than the size of the dot, the classical 
dynamics of the electron inside the dot is determined by the scattering from the dot's 
boundary. Due to small irregularities in the dot's shape, the electron displays chaotic motion, 
and its quantum transport through the dot can be described by statistical S'-matrix theory.0 
Since the Coulomb blockade regime is dominated by resonances, the conductance peaks can 
be used to probe the chaoticity of the underlying resonance wavef unctions. A statistical 
theory of the conductance peaks was originally developed in Ref. ||. By using i?-matrix 
theory,!'0 the conductance peak amplitude was expressed in terms of the electronic resonance 
wavefunction across the contact region between the dot and the leads. When the dynamics of 
the electron inside the dot is chaotic, the fluctuations of the wavefunction inside the dot are 
assumed to be well described by random matrix theory (RMT). In Ref. [8] the conductance 
distribution was derived in closed form for one-channel leads. These results were rederived 
in Ref. |IT|.and later extended to the case of two-channel leads in the absence of time-reversal 
symmetrytl through the use of the supersymmetry technique. Ill However, the calculations 
required by this technique become too complicated to apply in the general case of any 
number of possibly correlated and/or non-equivalent channels. 

The conductance distributions for one-channel leads were recently measure and 
found to be in agreement with theory for both cases of conserved and broken time-reversal 
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symmetry. This indicates that the dephasing effect, which plays an important role in open 
dots,@0 is of little importance for closed dots. 

In this paper we discuss in detail the width and conductance peak distributions for 
leads with any number of channels that are in general correlated and non-equivalent. Exact 
closed expressions for these distributions are derived for both cases of conserved and broken 
time-reversal symmetrylll We find that these distributions are entirely characterized by 
the eigenvalues of the channel correlation matrices M l and M r in the left and right leads, 
respectively. The strength of our approach is in its simplicity, since it relies solely on 
standard RMT techniques. To test our predictions we compare our analytical findings to 
numerical simulations of a chaotic dynamical model, the conformal billiard.lll Statistical 
width and conductance distributions of one-channels leads were recently studied in detail in 
this model.0 Although our paper deals mainly with ballistic dots whose classical dynamics 
is chaotic, our results should also be valid in the diffusive regime of weakly disordered dots, 
where random matrix theory is applicable. 

We note that under certain conditions the partial width is analogous to the wavefunction 
intensity at a given point. Therefore our width distributions can also be tested by microwave 
cavity experiments,0'il where the intensities are measured at several points that are spatially 
correlated. 

The plan of the paper is as follows: In Section II we briefly review the conductance in 
quantum dots in their Coulomb blockade regime. In Section III we discuss the statistical 
model and derive analytic results for the partial and total width distributions in each lead, 
for the channel correlation matrix and for the conductance distribution. We investigate the 
variation of these distributions as a function of the number of channels and their sensitivity 
to the degree of correlations between them. Those findings are compared in Section IV with 
numerical results obtained for the conformal billiard. Finally, in Section V we discuss the 
validity of our assumptions in the the context of typical experiments. 



II. CONDUCTANCE IN QUANTUM DOTS 

In this section we briefly review the formalism and introduce the notation used through- 
out this paper. In particular, we express the conductance peak heights in terms of the 
channel and resonance wavefunctions of the dot. 

For T <C kT <C A, which is typical of many experiments,! the observed on resonance 
conductance peak amplitude is given byffl 

G\ = — —— 9a with q x = , A A , (1) 
h 2kT yx yx T l x + ' K J 

where F x is the partial decay width of the resonance A into the left (right) lead. Since each 
lead can support several open channels we have = ^2 c T l } x , where T 1 ^ is the partial 
width to decay into channel c in the left (right) lead. 

In the .R-matrix formalism,^ the partial widths are related to the resonance wavefunction 
inside the dot. More specifically, introducing the partial amplitudes 7 c a, such that r cA = 
|7 c a| 2 , one can write 
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7cA = \F^ fdS<S>* c (r)V x (r). (2) 

Here ^\(r) is the A-th resonance wavefunction in the dot, $ C (V) is the transverse wavefunc- 
tion in the lead that corresponds to an open channel c, and the integral is taken over the 
contact area between the lead and the dot. k c and P c are the longitudinal wavenumber and 
penetration factor in channel c, respectively. 

Eq. @ shows that the contributions to the partial width amplitude from the internal 
and external regions of the dot factorize. The information from the region external to the 
dot is contained in k c and P c . These quantities are determined by the wave dynamics in 
the leads and are non-universal. They affect the average widths and enter explicitly in the 
correlation matrix M. However, the fluctuation properties of the conductance are generic 
and depend only on the statistical properties of the electronic wavefunction at the dot-lead 
boundary inside the barrier region. 

A different physical modeling of a quantum dot assumes point-like contacts and each 
lead is composed of several such point contacts. □0 In this model the conductance is also 
given by_(P) with each point contact r c considered as one channel. The corresponding partial 
width isB 



7cA = \P^*A(r c ) , (3) 

where A is the area of the dot, A is the mean spacing and a c is a dot-lead coupling parameter. 

Both models can be treated by our formalism. This becomes apparent after the following 
considerations. A resonance eigenfunction with eigenenergy E = E\ can be approximated 
by an expansion in a fixed basis p M of wavefunctions with the given energy E inside the dot 

* A (r) = 5> v p,(r). (4) 

The sum over /i is truncated at iV basis states, where N is large and determined by precision 
requirements. The partial width in channel c can then be expressed by the scalar product 

lc\ = (<t>c\ll>\) = <^/>A M , (5) 



where 

for the extended leads model, and 



l h 2 k c p c 



m 



dS$* c (r)p,(r) (6) 



^y^pyrc) (7) 

for the point contact model. Thus, we are led to similar formulations of both the extended 
leads and point-like contacts problems; in the corresponding iV-dimensional space the partial 
width amplitudes of a level are simply the projections of its corresponding eigenstate vector 
?/>a on the channel vectors <f> c . The only difference between the two models is the explicit 
expression for the channel vector C . We note that the scalar product (||) (that will be used 
throughout this paper) is different from the original scalar product defined in the spatial 
region extended by the dot. 
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III. STATISTICAL MODEL 



Due to the irregularity of the dot's shape, the motion of the electron inside the dot is 
expected to be chaotic. In Ref. |8| we have developed a statistical theory of the conductance 
peaks by assuming that the vectors ip\ that correspond to the resonance wavefunctions inside 
the dot have the same statistical properties as the eigenvectors of a random matrix ensemble. 
Here we study the limits of conserved time-reversal symmetry, corresponding to the Gaussian 
Orthogonal Ensemble (GOE), and of broken time- reversal symmetry, corresponding to the 
Gaussian Unitary Ensemble (GUE). The transition from one symmetry to another occurs 
when an external magnetic field is applied. The width distribution (or equivalently the 
wavefunction intensity distribution) was derived in the crossover regime between symmetries 
for the case of one channel leads only.il 



A. The Joint Distribution of Partial Width Amplitudes 



In RMT the eigenvector if) = (ipi,ip2, ■ ■^^n) (here and in the following we omit the 
eigenvector label A) is distributed randomly" on a sphere P{tf>) oc <5(EjU I^P ~ !)• The 
joint distribution of the partial width amplitudes 7 = (71, 72, • • • , 7a) for A channels is then 
given by 



P{l) 



V{/3N/2) 



7T 



(3N/2 



Dm 



fp(7e-<0eW) 



.c=l 



N 



(8) 



yit=l 



where D[if>] = nJU #a» for the GOE and D[if>] = nj=i #£<%/ 2 ™ for the GUE. To 
evaluate @ we transform the A channels to a new set of orthonormal channels C 



4>c = E 4>d F dc with (0J0 C /) = S c 



(9) 



We then take advantage of the invariance of the corresponding Gaussian ensemble under an 
orthogonal (unitary) transformation to rotate the eigenvector if} such that its first A com- 
ponents are along the new orthonormal channels. Denoting by O the orthogonal (unitary) 
matrix whose first A rows are the orthonormal vectors c (c = 1, . . . , A), we change variables 
in (§) to = Y.u CVVv Using ^ c = (0 C |V) we find 



P{l) 



T(PN/2) 



71 



/3N/2 



N 



At=A+l 



Vc=l 



f[*(7c 

Lc=l 



N 



X 



s EhA c | 2 + E \A 



,c=l 



M=A+1 

The integration over these first A components is now easily done and gives 

T(PN/2) 



P{l) 



W 2 I det F\ 



//J[^]5(7 t 7+ E l^| 2 -l] , 



(10) 



(ii) 
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where % = (<fi c \il>) = J2c> Ic'Fj^* are the partial widths to decay to the new channels and the 
metric is as before but excluding the first A components of if). Finally, the latter integral is 
easily done by introducing spherical coordinates in the N — K dimensional space. We obtain 

n~f) ~ ^A/ap {m _ A)/2) |detF| [1 7 * ) 7j • W 

For A C and in the limit A" — > oo, we recover a simplified expression 

P( 7 ) = (det M) _/3 / 2 e -i^ A/ "^ , (13) 
where the matrix M = (A^F^F) -1 is just the metric defined by the original channels 

M cd = ^{<t> c \<t>c<) ■ (14) 

The distribution (fi3| ) is normalized with the measure D[i\ = nli ^7c/2tt for the GOE 
and D[y] = Ilc=i d , y*d'y c /2ni for the GUE. Note that for both ensembles the joint partial 
width amplitudes distribution is Gaussian, the main difference being that the partial ampli- 
tudes are real for the GOE and complex for the GUE. Such a Gaussian distribution is also 
obtained by assuming that the distribution is form-invariant under an orthogonal (unitary) 
t r ansf ormation .0 

It follows from ( |13|) that the matrix M is just the correlation matrix of the partial widths 

M cd = t!t7. (15) 

In general the channels are correlated and non-equivalent, i.e. non-equal average partial 
widths. According to (fb3| ) this is equivalent to assuming channels that are non-orthogonal 
and have non-equal norms. 



B. The Channel Correlation Matrix M 



We shall now derive explicit expressions for the correlation matrix M in a chaotic quan- 
tum dot. Using Eq. fll4|) and the definition of the scalar product fl5|) we find 




We first discuss the case where there is no magnetic field so that the motion inside the 
dot is that of a free particle. Therefore, a resonance eigenstate inside the dot at energy 
E = h 2 k 2 /2m can be expanded in a basis of free particle states at the given energy E. Since 
RMT is applicable on a local energy scale, this is the fixed basis p M for which the eigenvector 
coefficients ip^ are distributed randomly (on the sphere). Using polar coordinates, such a 
basis of free waves is given by oc J^(kr) exp(z/i6 l ) with \x = 0, ±1, ±2, . . ., where J M are 

Bessel functions of the first kind. Denoting by A^ the number of such waves on the energy 
shell, we find 
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^£P»P>') = iE«W^ M = ^o(*|r-r'|) , (17) 

where we have used the addition theorem for the Bessel functions.il A similar relation holds 
if we choose a plane waves basis p M (f) = A^ 1 ^ 2 exp(ik^ ■ r) at a fixed energy h 2 k 2 /2m 
but with random orientation of fc M and use the integral representation of Jo- With help of 
Eq. ( |I?D we obtain for the correlation matrix 

M cc , = A-^Jk c k d P c P d Jdsj dS' <5>* c (r)J (k\r - r'|)$ c (r') , (18) 

for extended leads, while for the point contact model we find 

n A 

M cc , = ^—J {k\r c -r' c \). (19) 

7T 



We remark that Eq. (|I9[) is equivalent to C(k\Ar\) = ^*(r)^(r')/\^(r)\ 2 = J (k\r - r'\). 
This result was first derived in Ref. ^9] based on the assumption that the Wigner function 
of a classically chaotic system is microcanonical on the energy surface, and recently studied 
extensively in the Africa billiard.0 However, in these references the average is taken for a 
fixed wavefunction over a local region around (r + r')/2. 

When an external magnetic field B is present, the electronic classical underlying dynam- 
ics undergoes a transition from chaotic to integrable as the field gets stronger, regardless 
the shape of the billiard. In this paper, however, we only discuss the case of weak fields 
for which the motion is chaotic, and we are interested in the transition from orthogonal 
to unitary symmetry. While in the unitary case the wavefunctions become complex, the 
arguments that lead to Eq. (pj]) are still valid and the wavefunction correlator C(k\Ar\) is 
unchanged. 

The wavefunction correlation C(k\Ar\Yhas been also derived for weakly disordered sys- 
tems using the supersymmetry techniquecH in the unitary and orthogonal symmetries. In 
addition Ref. [H| derives the joint probability distribution for the intensity of an eigenfunction 
at two different points. We remark that the joint distribution of the wavefunction amplitude 



at A points r c is a special case of (|T3| ) obtained for 7 cA = ^\(r c ) (see the point contact case 
(0) except that the points r c can be chosen anywhere within the dot and not only on the 
boundary). We then obtain 



P(* A (n),* A (r 2 ),...,* A (r A )) = (detM)"^ 2 exp 



cc'=l 



, (20) 



where M cc > = A 1 J (k\r c — r' c \). The distributions of Ref. [31] are then easily obtained from 
(0) when A = 2.11 



C. Total Width Distribution 

We calculate next the total width distribution P(T) in a given lead that supports A 
channels and is characterized by a correlation matrix M. Although this quantity is not 
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directly measurable in experiments with quantum dots, it appears very often in resonant 
scattering by complex objects. We remark that for a dot with reflection symmetry T = 
r r = T the conductance peak g in ([I]) is proportional to T. Such dots are, however, difficult 
to fabricate. 

Using T = J2 C 1 7c 1 2 = 7 7, the characteristic function of P(T) is given by 

/oo /*oo 
dTexp(itT)P(T) = D [7] exp(i<7 + 7) P(j) . (21) 
-00 J —oo 

Since P("y) is a Gaussian, we readily obtain P(t) = [det(J — 2iMt/(3)]~ 13 ^ 2 . The distribution 
itself is then given by an inverse Fourier transform 

1 /-oo p—itr 

P(T) = — dt ^ . (22) 

2tt7-oc [det(/ - 2ttM/(3)f /2 



The matrix M is Hermitean and positive definite (since x'Mx = \x ■ 7I 2 > for any 
x 7^ 0) and therefore its eigenvalues w 2 are all positive. According to (P%D, P(T) depends 
only on w 2 . This is a consequence of the invariance of T under a orthogonal (unitary) 
transformation of the A partial width amplitudes. 

We first discuss the simpler GUE case, for which the integrand has poles —i/w 2 along the 
negative imaginary axis. Taking a contour integration along the real line and a half-circle 
that encloses all the poles in the lower half of the plane, we can calculate ( p2[) by residues. 
Assuming that all eigenvalues of M are non-degenerate, the poles are all simple and we find 



GUE 



(T) 




n ( 2 ,..2) 

c'+c W d W c 



e~ r/w " . (23) 



The distribution Pgue{F) given by (^) must be positive, which can be directly verified by 
using the concavity of the exponential function. 

For two channels (A = 2) which are in general non-equivalent (M n 7^ M22) and correlated 

(Mi 2 0), the eigenvalues are given by w 2 l2 = (M 11 + M 22 )/2±J ((M u - M 22 )/2) 2 + \M 12 \ 2 



Then, Eq. ( |23[ ) reduces to 



P a= 2 (£) - 2a+ e -*$m-\m sinh [ 2a+ V a - + l/|2 f) (24) 
P GUE {V) - ^—-—e smh^ x _^ 2 FJ , (24) 

where T = T/T is the width in units of its average value, / = M 12 / y /M 11 M 22 measures the 
degree of correlation between the two channels and a± = 1/2 ^M u /M 22 ± \J M 22 /M u ^ are 
dimensionless parameters such that for equivalent channels a + — 1 and a_ = 0. In the latter 
case, we reproduce the result of Ref. [12]. 

For degenerate eigenvalues, we can calculate (p2|) by using the residue formula for 
higher order poles. Alternatively we can slightly break the degeneracy of the eigenval- 
ues by 77 and take the limit 77 — >• 0. For example, for two channels Eq. ( ^3|) gives 
-P(r) = ^e~ r / w i — e~ v /' w 2^ /{w 2 — wf). By taking w 2 = wj + 77, in the limit 77 — * we 
recover 
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(25) 



which is the x 2 distribution in four degrees of freedom. More generally, when all A channels 
are uncorrelated and equivalent (M = w 2 I) we recover the well-known \ 2 distribution in 2A 
degrees of freedom! 

p ^ (^) = ^^^^)! ^A " e " / " , • (26) 

We have denoted this limiting distribution in (|26|) by Pque as ^ wm serve as our reference 
distribution against which to compare the distributions in the general case of correlated 
and/or inequivalent channels. 

For the GOE case, the integral of Eq.(|22|) is more difficult to evaluate since the sin- 
gularities of the integrand along the negative imaginary axis t = —ir are of the type 
(r — l/2w 2 )~ 1 / 2 . In this case the semi-circle part of the contour (in the lower half of 
the plane) is deformed to go up and then down along the negative imaginary axis so as 
to exclude all the singularities. When going around a singularity of the above type the 
function changes sign. Therefore, after sorting the inverse eigenvalues of M in ascending 
order wf 2 < w^ 2 < . . ., we have 

1 /„ 1 \ „ fV^im e~ FT , s 

wn = ^ n- eL 2 dr /n ^-u T i \ tta h ;t (27) 

where for an odd number of channels A, we define l/2w\ +l — ► oo. The integrand of each 



term on the r.h.s. of fl27|) is singular at the two endpoints of the integration interval, but 
this singularity is integrable. For the case of two channels that are in general non-equivalent 
but correlated, Eq. (|27|) reduces to 



v/i-l/l 2 \ l ~W 



where / and a± are defined as before (see following Eq. ([24])) and I is the Bessel function 
of order zero. The case of equivalent channels is obtained in (|28|) by substituting a + = 1 
and a_ = 0. 

The reference distribution Pgoei defined as before for the case where all A channels 
are equivalent and uncorrelated, is found directly from (|22|) to be the \ 2 distribution in A 
degrees of freedom!!'! 

^' w^-Dr 1 '^' (29) 

The top panels (a and b) in Fig. show the width distributions for a two-channels lead in 
the GOE statistics. The left panel is for equivalent channels (M22/M11 = 1) and for various 
degrees of correlations / = 0.25, 0.5, 0.75, and 0.95. The right panel is for uncorrelated 
(/ = 0) but non-equivalent channels: M^jMw = 2,3,4, and 5. The bottom panels (c and 



9 



d) in Fig. [I] are similar to (a and b) except that they correspond to the GUE case. We note 
that in all figures we display as T the normalized total width T/T. 

The correlation matrix in the point contact model is fully determined by fc|Ar| and 
the number of channels A. The left panels in Fig. H show the GOE width distributions 
for fc|Ar| = 0.25,1,4 and for different number of channels A = 2,4, and 6. The right 
panels of Fig. ^| show similar results but for the GUE statistics. The deviation of the width 
distribution from the reference distribution P'°'(r) which corresponds to equivalent and 
uncorrelated channels (dashed lines in Fig. 0) becomes larger as the number of channels 
increases for a given fc|Ar|. 



D. Conductance Peaks Distribution 



To calculate the conductance distribution P(g) in the general case, we assume that the 
left and right leads are far enough from each other and thus uncorrelated. @ The left and 
right leads are characterized by their own correlation matrix M l and M r , respectively. Under 
this assumption 



P(g) = JdT l dT r 5(g 



Ipr 



rT 



p(r*)P(r r ) 



where P{T) is given by ( |23| ) in the unitary case and by fl2~T|) in the orthogonal case. 
The distribution P(g) can be evaluated by the following identity 

/■oo roc / 

J dT 1 J dT 2 e- r '/ 6l e- r ^ S(g - 




(30) 



(31) 



provided 61,62 > 0. To obtain this identity we have used the integral representation of the 
Bessel function K v {z) = l/2(z/2) v J^dtt-^e^ 2 ^ . 
For the unitary case, Eqs. ( |2~3"D and (|3lD give 



n ( 2 .,2) n ( 2 ...2) 

cVc V d V c d'^d W d> W d 



xe ? ™5 



Kr 



2g_ 

V c W d 



2g 



2 \w d v c ) \v c w d 



(32) 



where v\ and w\ are the eigenvalues of the left and right lead correlation matrices M l and 
M r , respectively. The previous published resultsi§ are special cases of Eq. (|T5|) for one 
channel leads with r = T (i.e. v% = w{), while the distribution of Ref. |1^ is obtained for 
two (equivalent) channels leads whose matrices are related by an overall asymmetry factor 
M r = aM l . 

A similar calculation for the orthogonal limit gives 
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Fig. |3| shows the GOE (left) and GUE (right) conductance peak distribution (^) and 
(|32|), respectively, for symmetric A-point leads with fc|Ar| = 0.25, 1,4 and for A = 2,4, and 
6 (the same cases shown in Fig. §). I n analogy to P(T), all figures depicting P(g) display 
the normalized conductance g defined as g/g. By comparing Fig. |3]with Fig. ^|we conclude 
that, as A increases, the conductance distribution shows stronger deviation from its limiting 
case of uncorrelated equivalent channels (dashed lines) than the width distribution does. 

Fig. H shows the case of asymmetric leads for the asymmetry factor a = 1 and 10, 
for four-point leads with fc|Ai*| = 1 and for the orthogonal (a) and the unitary (b) limits. 
P(g) is not very sensitive to the leads asymmetry and a large value of a is needed to see 
significant variation from the symmetric leads case. In the limit a — > oo, one can neglect 
the smaller width in ([!]) and the conductance peak g is proportional to the partial width 
in the dominating lead. In this limit P(g) is reduced to P(T) shown by the dashed lines in 
Fig. |j). The asymmetry effect becomes larger for an increasing number of channels. This 
effect can be noticed by comparing the GOE and GUE cases, since for the same number of 
physical channels A the GUE has a larger number of "effective" channels. 



IV. DYNAMICAL MODEL 

To test the RMT predictions for the statistical distributions, we modeled a quantum dot 
by a system whose classical dynamics is chaotic. The model is the conformal billiardjllll 
whose shape is defined by the image of the unit circle in the complex z-plane under the 
conformal mapping 

. , z + bz 2 + ce i5 z 3 . . 

w(z) = . . (34) 

V ' VI + 2b 2 + 3c 2 V ; 

The parameters b, c and 5 control the billiard shape. Eq. (|34|) ensures that area A enclosed 
by w(z) is normalized to it and is independent of the shape. We analyze the case b = 0.2, 
c = 0.2 and 5 = ir/2, for which the classical phase space is known to be chaotic.129 We have 
verified that the corresponding spectrum exhibits GOE-like spectral fluctuations (we used 
300 converged levels by diagonalizing a matrix of order 1000). This is demonstrated in Fig. 
[5], where the nearest-neighbors level spacing distribution P(s) and the A3 statistics, which 
measures the spectral rigidity, are shown. 

To investigate the effect of an external magnetic field, we consider the same billiard 
threaded by an Aharonov-Bohm flux line,iH which does not affect the classical dynamics. 
The flux is parametrized by $ = a&o where $0 is the unit flux. We use the same set of 
values for b, c, and 5 to insure classical chaotic motion, and choose a = 1/4 for maximal 
time-reversal symmetry breaking. The statistical tests shown in Fig. [5] confirm that this 
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choice of a corresponds to the unitary limit . We remark that the A3-statistics is a better 
measure to distinguish between the GOE and GUE cases than the level spacing distribution 
P{s) (used in Ref. |20|.) 



A. Spatial Correlations 



The eigenfunction amplitude correlation C(k\Ar\) = \E'*(r)^(r / )/|\l'(r)| 2 was recently 
investigated thoroughly for the conformal billiard.^ The results agree fairly well with the 
theoretical prediction, namely C{k\Ar\) = J^{k\r — r'\)^ if one averages over the orientation 
of Ar. This result is obtained based on semiclassical arguments and the eigenfunctions stud- 
ied in Ref. ^0| were chosen accordingly to be highly excited states (deep in the semiclassical 
region) . 

In order to apply this result to quantum dots, further considerations are in order. First, 
a typical semiconductor quantum dot in the submicrometer range contains several hundred 
electrons, and it is therefore not obvious that the eigenstates around the Fermi level are 
necessarily semiclassical. Second, scars associated with isolated periodic orbits give correc- 
tions to C(k\Ar\) which depend on the orientation of Ar and are of order 0(fr 1//2 ). The 
fluctuations of the spatial correlation of the billiard eigenfunctions were recently studied^ 
and found also to be suppressed by 0{h 1 ^ 2 ). These corrections are negligible if one averages 
over all orientations around a given point r, keeping the modulus |Ar| fixed, but this is 
difficult to implement experimentally. At a fixed orientation the fluctuations of the spatial 



correlations seem to be rather small if fc|Ar*| 3 so that (|T9|) is a good approximation. For 
larger values of fc|Ar|, there could be significant fluctuations from ([19]) but in this region 
the width and conductance distributions are closer to their limiting case of independent 
channels and are not very sensitive to the exact correlations. 

Our results were obtained by using the billiard eigenfunctions with Neumann boundary 
conditions where the normal derivative of the wavefunctions vanishes on the boundaries. 
We analyze eigenfunctions in the vicinity of the 100 t/l excited level which resembles the 
experimental situation. By moving the points around the circle we generate more statistics 
and average over orientations. The results are shown in Fig. |6| where the correlations in 
the model (solid line) compare well with the theoretical result (dashed line) for both cases 
with and without magnetic flux. The agreement is fair, particularly for /c|Ar| < 5. For 
fc|Ar| ^> 1, the deviations from the theoretical value of C(k\Ar\) are not important since 
the channels are weakly correlated and the distributions are very close to those describing 
uncorrelated channels. Thus, corrections to our analytical findings should not be large, as 
is supported by the numerical evidence presented below. 

In our model studies we imposed Neumann boundary conditions around the entire billiard 
and not just at the dot-lead boundary. To mimic the experimental situation we would 
have to impose mixed boundary conditions,^ which makes the calculations much more 
computationally intensive. However, we now argue that our simplified situation still provides 
reasonable results. For extended leads, the length of the dot-lead contact region D must 
satisfy kD ^> 1 in order to support open transverse channels (in dots containing several 
hundred electrons). Therefore, deviations from C(k\Ar\) at the edge of the dot-lead contact 
region (where our boundary conditions are unrealistic) are averaged out. For point-like 
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contacts, the physical picture is that the conductance is probing the wavefunction in the 
vicinity of the constriction (the region that couples the dot to the external lead). We then 
need to know the characteristic properties of the wavefunction inside the dot where our 
model is quite satisfactory 

B. Coupling to Leads and Distributions 

We first studied the point-like contacts model by describing the lead as a sequence of A 
equally spaced points on the boundary of the billiard (in the to-plane) . According to Eq. ( |19| ) 
the correlation matrix M is then completely determined by fc|Ar| ~ kS9\w'(r = 1,9)\ (where 
59 is the angle that spans the arc between two neighboring points in the z-plane) and A. 
In this model it is easy to generate strongcorrelations by choosing the points close enough, 
unlike the (discretized) Anderson modelll3 where the channels are weakly correlated even 
if the lead is composed of nearest neighboring points. The eigenvalues w\ are found by 
diagonalizing the matrix M. 

In Figs. [7] and || we compare for the unitary and orthogonal limits, respectively, the total 
width distribution P(T) obtained by solving the conformal billiard (histograms) with the 
theoretical predictions (solid lines), for several values of fc|Ar| = 0.5, 1, 2 and A = 2, 4, 6. The 
distributions P^(T) for equivalent and uncorrelated channels are indicated by the dashed 
lines, and are just the x 2 distributions in A (2A) degrees of freedom for the GOE (GUE). The 
agreement between the model and the analytic RMT predictions confirms the validity of the 
statistical model for a chaotic dot. We observe from Figs. [7] and |8| that for the larger values 
of fc|Ar|, the distributions get closer to those for uncorrelated channels. This is consistent 
with the decrease in spatial correlations (see Fig. [7]). Another interesting observation is 
that, for a constant fc|Ar| (i.e. fixed correlations), the deviation from the limiting case of 
independent channels becomes larger with an increasing number of channels. 

Figs. |9] and [U] show a comparison between the theoretical conductance peaks distribu- 
tions for symmetric leads, as given by Eqs. fl3"2] ) and for the unitary and orthogonal 
cases, respectively, and those calculated for the Africa billiard with symmetric A-point leads 
(A = 2,4, and 6) and for different values of fc|Ar|. The dashed lines are again the limiting 
case of uncorrelated and equivalent leads. Observations that are similar to the ones made 
above for the width distributions, can be made with respect to the conductance peaks dis- 
tributions. Comparing the width and conductance peaks distributions, we note that the 
conductance distribution shows stronger deviation from its limit for uncorrelated equivalent 
channels than does the width distribution. 

We also studied extended leads by taking the contact region of the lead and the dot 
to have a finite length D rs \w'\A9 on the dot's boundary (in the w-plane) where w' is 
evaluated at the corresponding angle where the lead is located. In this case the channels 
are defined by the allowed quantized transverse momenta k c = irn c /D with n c — 1, 2, . . . , A, 
where A = mt[kD/7r]. To calculate the partial amplitude for the conformal billiard, the 
integral in Eq. @ (defined in the w-plane) is mapped into an integral along an arc in the 
z-plane which is spanned by an angle A9 

7cA = V2^ Le M IW ' {T = 1,9)1 ^(^(r = 1,9), (35) 
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where $ c (#) = y2/Dsm(K c \w'\9) are the transverse channel wavefunction and for simplicity 
we have set k c P c = 1. The resonance eigenfunction ty\ is given in terms of its expansion in 
e im6 (withm = 0,±1,±2,...) 



*x(r = 1,6) = A/" A £ 1 . = ' = ^ , (36) 

j y/nij-\tj-<xn 

where Af\ is a normalization constant, jj are the zeros of JL._ a \ and Cj are expansion 



coefficients as in Ref. 20 



To guarantee that the correlation matrix M in ( |lTj|) is the same for eigenf unctions of the 
billiard which belong to different energies, we choose D such that kD = constant and scale 
the partial amplitude (|j) by k 1 ^ 2 . The resulting matrix is 

kM cd = ( dO ( dd'\w\r = 1,0)| \w'(r = 1,9') \ 

2m kD JA6 J Ad 

xsm(^\w'\ey o (\w'\\e-6'\) S m(^\wY) ■ (37) 

This scaling is desirable in order to be consistent with the theoretical approach presented 
above, but experimentally it is vervhard to accomplish. Fortunately, this scaling of D is 
insignificant for present experiment S00 that deal with dots containing several hundred elec- 
trons Af. Indeed, from the Weyl formula we have k F oc A/" 1//2 so that Skp/k F = 5N/2N' 1. 
The latter inequality is obtained when we estimate 5AF to be the number of observed 
Coulomb blockade peaks (since each Coulomb blockade peak corresponds to the addition of 
one electron into the dot). The relative variation of kp is thus small and can be neglected. 

We find that the channels in the extended leads model are weakly correlated and that the 
average partial widths in the various channels exhibit a moderate variation. In such a case 
the total width distribution is not very different from the case of uncorrelated equivalent 
channels. Our model calculations for extended leads are shown in Fig. |Tl| and are in agree- 
ment with the RMT predictions for uncorrelated channels (dashed lines). An interesting 
effect is that with an increasing number of channels even small deviations in P(T) give rise 



to relatively large deviations in P(g), as can be seen in Fig. lid 



V. CONNECTION TO EXPERIMENTS AND CONCLUSIONS 

We have discussed both the cases of orthogonal and unitary symmetries. To relate 
to actual experimental situations, it is important to estimate the minimal strength of the 
magnetic field B c which ensures complete time-reversal symmetry breaking. For a ballistic 
electron9§@'H B C A oc y / r c , r /r$ , where r cr and r are, respectively, the time it takes the 
electron to cross the dot and the Heisenberg time r = h/A. For an electron at the Fermi 
energy B C A oc A/" -1 / 4 ^) where M is the number of electrons in the dot. The proportionality 
factor is non-universal and depends on the exact geometry of the dot. In a semi classical 
analysist^l'EJ it can be expressed in terms of classical quantities. For the dots used in some 
recent experiment s,0E! B c is of order of a few mT. Such small values of B c do not alter 
significantly the classical dynamics of the electron,0 and our assumption that the correlation 
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C(k\Ar\) is unchanged is justified. Nevertheless, these small variations in the magnetic field 
have appreciable quantum mechanical effects, i.e. the crossover from orthogonal to unitary 
symmetry. 

In conclusion, we have derived closed expressions for the width and conductance peak 
distributions in quantum dots in the Coulomb blockade regime. The main assumption is that 
the electron's dynamics is chaotic for a ballistic dot or weakly diffusive for a disordered dot. 
For given correlation matrices that characterize the left and right leads, these distributions 
are universal and distinct for conserved and broken time-reversal symmetry. While recent 
experiments have measured the conductance distributions in symmetric one-channel leads, 
it would be interesting to measure and compare with theory the conductance distributions 
in more general situations of dots with multi-channel leads. 
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FIGURES 



FIG. 1. Total width distributions P(T) for a two-channel lead. Panels (a), (b) correspond to 
the orthogonal symmetry and (c), (d) to the unitary symmetry. P(T) for equivalent but corre- 
lated channels with / = 0.25,0.5,0.75 and 0.95 are shown in (a) and (c), while uncorrelated but 
non-equivalent channels M22/-M11 = 2,3,4,5 are shown in (b) and (d). Descending values of / (or 
M22/M11) correspond to distributions that extend to larger values of F 

FIG. 2. Total width distributions for a A-point lead and fc|Ar| = 0.25, 1 and 4 in a quantum 
dot with orthogonal symmetry (GOE) and with unitary symmetry (GUE). (a) GOE; A = 2; (c) 
GOE; A = 4; (e) GOE; A = 6; (b) GUE; A = 2; (d) GUE; A = 4; (f) GUE; A = 6. The dashed 
lines correspond to uncorrelated and equivalent channels. Increasing values of fc|Ar| correspond 
to curves which approach the case of uncorrelated channels. 

FIG. 3. Same as Fig. |2|but for the conductance peak distributions P(g) in dots with symmetric 
A-point leads. 

FIG. 4. Conductance peak distributions P(g) for asymmetric four-point leads with fc|Ar| = 1 
and an asymmetry factor of a = 1 and a = 10. (a) GOE; (b) GUE. The dashed curves describe 
the limit a — > 00 where P(g) reduces to P(T). 

FIG. 5. The nearest neighbors level spacing distribution P(s) and the A3-statistics for the 
conformal billiard with b = 0.2, c = 0.2 and 5 = ir/2. We consider the states between the 50th and 
the 350th. Left: no magnetic flux ( a = 0). Right: with magnetic flux of a = 1/4. 

FIG. 6. The spatial wavefunction correlation C(k\Ar\) calculated for the conformal billiard 
(squares) compared with the theoretical prediction Jo(A;|Ar|) (solid line). Panel (a) displays the 
case where a = and (b) corresponds to a = 1/4. 

FIG. 7. Total width distributions P(T) for the unitary case for several values of 
k\ Ar| = 0.5, 1, 2 and for various number of channels A = 2, 4, 6. The solid lines are the theoretical 
distributions (^), while the dashed lines correspond to uncorrelated and equivalent channels. The 
histograms are the results from the Africa billiard (6 = 0.2, c = 0.2,5 = 7r/2) where a magnetic 
flux line (a = 1/4) breaks the time-reversal symmetry. 

FIG. 8. Total width distributions P(T) for conserved time-reversal symmetry (GOE). Con- 
ventions are as in Fig. ^, with solid lines describing the orthogonal prediction (27). 



FIG. 9. Conductance peak distributions P(g) for the unitary symmetry. The histograms 
display the results obtained from the Africa billiard for symmetric A-point leads and different 
values of fc|Ar|. The solid lines are the RMT prediction (|32|) and the dashed lines correspond to 
uncorrelated and equivalent channels. The cases presented are the same as in Fig. ^. 

FIG. 10. Same as Fig. |9| but for the orthogonal symmetry where the theoretical distribution 
is given by Eq. (p3|) . 



18 



FIG. 11. Comparison of conformal billiard results for extended leads with mt[kD/ir] = 6 
(histograms) and theoretical predictions. The panels represent total width distribution P(T) for 
orthogonal (a) and unitary (b) limits and conductance distribution P{g) for orthogonal (c) and 
unitary (d) limits. The dashed lines correspond to uncorrelated equivalent channels. 
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